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Abstract 

We consider effects of Coulomb interaction in a granular normal metal at 
not very low temperatures suppressing weak localization effects. In this limit 
calculations with the initial electron Hamiltonian are reduced to integrations 
over a phase variable with an effective action, which can be considered as a 
bosonization for the granular metal. Conditions of the applicability of the 
effective action are considered in detail and importance of winding numbers 
for the phase variables is emphasized. Explicit calculations are carried out 
for the conductivity and the tunneling density of states in the limits of large 
g ^ 1 and small g <^ 1 tunnelling conductances. It is demonstrated for 
any dimension of the array of the grains that at small g the conductivity 
and the tunnelling density of states decay with temperature exponentially. 
At large g the conductivity also decays with decreasing the temparature and 
its temperature dependence is logarithmic independent of dimensionality and 
presence of a magnetic field. The tunnelling density of states for g ^ 1 is 
anomalous in any dimension but the anomaly is stronger than logarithmic in 
low dimensions and is similar to that for disordered systems. The formulae 
derived are compared with existing experiments. The logarithmic behavior 
of the conductivity at large g obtained in our model can explain numerous 
experiments on systems with a granular structure including some high Tc 
materials. 

PACS numbers: 73.23.Hk, 73.22.Lp, 71.30.+h 
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I. INTRODUCTION 



Study of disordered systems attracts a lot of attention. Even if the electron-electron 
interaction is neglected the problem is not simple. Nevertheless, by now a rather complete 
understanding for non-interacting systems exists. It is well known that systems with one- 
dimensional geometry must be insulators for any week disorder. Two-dimensional systems 
are also believed to be localized, whereas there should occur the Anderson metal-insulator 
transition in three-dimensional disordered systems. This scenario has been formulated on the 
basis of scaling arguments long ago^ and has been checked by numerous methods including 
diagrammatic expansions^ and non-linear cr-model calculations^ (for a review, see e.g. Ref.^). 

At weak disorder, one obtains so called weak localization corrections which arc divergent 
in ID and 2D (in 2D they are logarithmic). The infrared divergency is cut off either by 
the sample size or external frequency, or by the inverse inelastic time. Within the scaling 
theory the concrete mechanism of this cutoff does not matter. Due to the existence of the 
cutoffs flow diagrams stop at some point when changing parameters like the sample size 
or temperature. So, within this picture by increasing the temperature one cannot expect 
anything except cutting the weak localization corrections. 

If an electron-electron interaction is added additional corrections discovered by Altshuler 
and Aronov^ become important. They are also logarithmically divergent in 2D but remain 
convergent in 3D. For the systems with the electron-electron interaction one can also derive 
a proper cr-model^ and demonstrate the renormalizability. Within such an approach the 
role played by the temperature is, as for the non-interacting models, to cut the diverging 
infrared corrections. After they are cut off the system is assumed to become a conventional 
normal metal. Although granular systems with the electron-electron interaction were not 
considered explicitly, it is generally expected that there should not be a qualitative difference 
between them and "just disordered" materials. This is definitely true for systems without 
the electron-electron interaction. 

Therefore, a recent experimental observation^ came as a surprise. In the experiment^ 
measurements of conductivity of granular Al — Ge thick films were performed and several 
interesting effects have been discovered. These films consisted of Al grains embedded in an 
amorphous Ge matrix. The size of an Al grain was about 120A and the grains were at low 
temperatures in a superconducting state. 

Destroying the superconducting pairing by a magnetic field up to 17T the authors could 
study, in particular, properties of the normal state. Some features of the normal state related 
to a negative magnetoresistance due to superconducting fluctuations have been discussed 
recently® but another unusual observation remained unexplained. 

The authors of Rcf.'' found in some samples a peculiar temperature dependence of the 
conductivity. Samples that had a high room temperature resistivity (this corresponds to a 
weak tunneling between the grains) showed an exponential decay of the conductivity as a 
function of temperature. This behavior is typical for insulators and has been interpreted in 
Rcf.^ in this way. Samples with larger intergranular couplings did not show any exponential 
decay but the resistivity did not saturate at low temperatures and the authors described its 
temperature dependence by a power law 

R = AT-'^, (1.1) 
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with a — 0.117. Apparently, with such a small value of a a logarithmic temperature depen- 
dence 



A{\- alnT) 



(1.2) 



(obtained by expansion of Eq. (1.1) in a) could describe the experimental data as well. 

Both the power law dependence, Eq. (1.1), and the logarithmic dependence, Eq. (1.2), 
were not expected from the theories of disordered metals. Although the logarithmic depen- 
dence is natural of ID films, it is relevant to emphasize that the array of the grains was 
three dimensional and the behavior described by Eqs. (1.2) was observed at high magnetic 
fields. This definitely excludes an interpretation of the logarithmic or power law behavior 
in terms of weak localization corrections typical for 2£)^'^'^. 

In this paper we try to explain the transition from the exponential temperature behavior 
at small g to the logarithmic dependence of Eq. (1.2) assuming that the temperature is not 
very low, which actually corresponds to the experiment^ where the lowest temperature was 
around O.SiiT. We consider a model of granular metals at not very low temperature and 
demonstrate that changing the dimensionless tunneling conductance g one can have either 
exponential temperature dependence of the resistivity at small < 1 or the logarithmic 
behavior, Eq. (1.2), at large ^ 1. It will be shown that the result is applicable for any 
dimensionality of the array of grains, which contrasts usual logarithmic corrections due to 
interference effects typical for ID only. 

Study of the regime of not very low temperatures of the granular metals has started 
recently in Ref.^ where it was demonstrated that the regime of the temperatures 



where b is the mean level spacing, is clearly different from the regime of lower temperatures. 
In this regime, well known weak localization effects are suppressed but, nevertheless, the 
system may exhibit a non-trivial behavior. 

Eq. (1.3) is written for g ^ \. In this limit, one can study the model by summation 
of diagrams containing impurity lines, as it has been suggested in Ref.^. Non-trivial effects 
arise from a renormalization of the Coulomb interaction by impurities. Actually, a non- 
trivial behavior originates from the diagrams of the Altshuler-Aronov type^ modified for the 
granular system. Therefore, for large (7 3> 1 the theory presented in this paper is equivalent 
to the summation of the Altshuler -Aronov diagrams. However we use an effective phase 
functional, which gives us an opportunity to consider not only the first order in the Coulomb 
interaction but to go into higher orders by writing renormalization group equations. 

If the tunneling conductance is not very large, fif < 1, the inequality (1.3) should be 
replaced by the inequality 



In the limit < 1, it is not possible to get any results by summation of conventional diagrams 
and one should use the phase functional. A very important result of Ref.^ is that in the 
limit of not very low temperatures, Eqs. (1.3,1.4), the phase functional has a form similar 
to that suggested by Ambegaokar, Eckern and Schon (AES) long ago^° for a description of 
quantum dissipation^^ in metals. 



T » gb, 



(1.3) 



(1.4) 
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At lower temperatures T < gS one must take into account low energy diffusion modes, 
and the system should be described by a more complicated non-linear cr-model containing 
both the phases and matrices Q. Of course, in the limit T — > the AES functional is 
no longer applicable and the notion of "quantum dissipation" looses its sense. In other 
words, the "quantum dissipation" corresponds to a limit of not very low temperatures for 
the system under consideration. 

It is relevant to mention that the model of a granular metal may roughly describe dis- 
ordered systems with a low electron concentration. One can imagine that in such systems 
electrons spend a long time in some kind of traps or puddles and tunnel between them with 
a small rate. Then, such a situation resembles the granular systems. 

A short account of the ideas of the present work has been presented in a recent paper^^ 

The article is organized as follows: 

In Sec. II we formulate the model and derive a free energy functional containing phases. 
We discuss how one should integrate over the phases accounting for winding numbers. In 
Sec. Ill we consider conductivity in the limit of large conductances g ^ 1 and obtain the 
temperature behavior, Eq. (1.2). For calculations we use both the perturbation theory and 
renormalization group techniques. In Sec. IV we consider the tunneling density of states 
in the limit g ^ 1. The final result for a two dimensional array of grains, Eq. (4.9), has 
the same form as the corresponding formula for a disordered "homogeneous" metal. In Sec. 
V we consider the limit of small g <^ 1 and show how one can carry out summation over 
winding numbers. Physical quantities like the conductivity and the tunneling density of 
states are shown to be exponentially small in temperature. In Sec. VI we discuss the results 
and make a comparison with the experiment. 



In the present work, we consider a simplified model where metal grains form a regular 
lattice. The grains are disordered due to impurities or irregular boundaries. Each grain is 
separated from its nearest neighbors by tunneling barriers. We assume that the main con- 
tribution to the macroscopic resistivity of the granular system comes from the intergranular 
tunneling. 

The Hamiltonian describing the model is chosen as 



II. CHOICE OF THE MODEL. PHASE FUNCTIONAL. 



H = Ho + Ht + H, 



Cl 



(2.1) 



where Hq is the one-electron Hamiltonian of isolated grains including disorder within the 
grains 




where U (r) is a random potential. 

The tunneling of the electrons between the grains is given by 




(2.2) 



i,j,a,a' 
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where the summation is performed over the states a. a' of each grain (spin is conserved) 
and over neighboring grains i and j. The possibihty to tunnel from the state a to an 
arbitrary state a! of other grains introduces an additional disorder resuhing in a finite 
tunnel conductance. 

The term He in Eq. (2.1) describes the charging energy 



is the excess number of electrons in the i-th grain, is the dimensionless local potential 
and Cij is the capacitance matrix. The integration over r in Eqs. (2.1, 2.4) is performed 
over the grain with the coordinate i of the center and includes summation over the spin. 

Eq. (2.3) describes the long range part of the Coulomb interaction in the limit of weak 
disorder inside the grains and has been used in many works (for a review, see^'^). In principle, 
one could consider also a superconducting or magnetic part of the interaction within the dot 
but we assume that the grains are in the normal state. The long range part of the Coulomb 
interaction, Eq. (2.3), describes the classical charging energy. This interaction can lead to 
the Coulomb blockade and to insulating macroscopic properties of the system of the grains. 

Calculations with the Hamiltonian Eqs. (2.1, 2.3), can be replaced in a standard way 
by computation of a functional integral over anticommuting fields (ri,r). 

Although the model described by Eqs. (2.1, 2.3) contains only the long range part of 
the Coulomb interaction, it is still very complicated, because at very low temperatures in- 
terference becomes very important and one has to consider an interplay of localization and 
interaction effects. One could do this either using diagrammatic expansions^ or writing 
a non-linear a-model^. Both methods allow to consider the limit of large tunneling con- 
ductances g and the results are strongly dependent on the dimensionality. However, the 
behavior, Eq. (1.1) or Eq. (1.2) was not predicted for 3D in any of these works. 

The model, Eqs. (2.1, 2.3), becomes simpler if the temperature T is not very low such 
that low energy diffusion modes are damped. As it was discussed in a recent publication^, 
the granular metal can be well described at temperatures T ^ gd^ where b is the mean level 
spacing in a single grain, by the Ambegaokar, Eckern and Schon (AES)^*^ functional of the 
free energy. If 51 < 1, this condition should be replaced by T ^ 5. The limit of not very low 
temperatures not only simplifies the consideration but is interesting on its own because it 
leads to an unusual behavior of physical quantities and is easily accessible experimentally. In 
particular, we will see that changing the tunneling conductance g one may have a crossover 
from the exponential temperature dependence of the resistivity to the logarithmic behavior. 



We calculate the conductivity o {lS) using the Kubo formula and making an analytical 
continuation from Matsubara frequencies = 2'KinT to real frequencies u?-^. Within this 
formalism the conductivity a [uS) can be written in the form 




(2.3) 



In Eq. (2.3), 




(2.4) 



Eq. (1.2). 
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la 



d-2 



(2.5) 



where a is a vector connecting the centers of neighboring grains i and i + a, a — |a| , and d 



is the dimensionahty of the array. 

The function Kg^i in Eq. (2.5) can be written as 



where 



and 



-R'aa' (ri, rj; t) = riaa' (r^, Tj] t) - e^(5aa'^ij5 (t) < H^^ (t) + h.c. > 



Haa' (ri, rj; r) = - < tJ^ (r;, r) 4' (rj, 0) > 



4l(^)=E^i+a,*,i+a(^)4',i(r) 
a,a' 



(2.6) 



(2.7) 



(2.8) 



The abbreviation h.c. means a Hermitian conjugation. The tunnehng current operator 
■^a(i"i,T) entering Eq. (2.7) takes the standard form 

Ja (ri) = ie (^H^, - h.c^ (2.9) 

The first term in Eq. (2.6) corresponds to a "paramagnetic" contribution whereas the 
second one is of a "diamagnetic" origin. 

In order to reduce the calculation of physical quantities to a computation of correlation 
functions with the AES action we decouple the interaction term, Eq. (2.3), by integration 
over an additional Vi (r) 



exp 



{r (ri,r)t^(ri,r) rfr; - iV) (r) dr 



^E/ dTV,{T)a^s{T)^ DV 



(2.10) 



and then, following Refs.^^'^, remove this field from Hq by the gauge transformation 

^ (ri,r) ^ e-^^'M^i (ri, r) , 0; (r) = (r) . (2.11) 

This is not a trivial procedure, because the new fields ipa {t) should obey, as before, the 
fermionic boundary condition 



V'(ri,T) = -V'(ri,T + /3), P^l/T 



(2.12) 



Let us forget for a while about the tunneling between the grains and explain how one can 
proceed for the Hamiltonian Hq of a single grain. 
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Due to the boundary condition, Eq. (2.12), the field (r) cannot be completely removed 
from Hq and this should be done approximately. The field Vj (r) can be represented as a 
sum of a static Vqi part and periodic function (r) 

V,{t) = Vo, + V,{t), ['Vi{r)dT = (2.13) 

The static part Vq can still be arbitrary and we rewrite it as 

Voi = 27rTki + Pi, - ttT < Pi < ttT (2.14) 

where h = 0,±1, ±2,.... 

If we neglected pi we would be able to remove Vj (r) from Hq. In this case we would use, 
instead of the phases (f (r) from Eq. (2.11), the phases (r) defined as 

0i (r) = 0i (r) + 27rTA;iT, (2.15) 

where —oo < 0i (r) < oo, 0i (0) = 0i {/3) . 

It is clear that making the gauge transformation, Eq. (2.11), with the phases 0i (t) 
instead of <^i(r), preserves the antiperiodicity of the ^ (r), Eq. (2.12). 

The variable pi cannot generally be neglected. This term is not important only in the 
limit T ^ S. In order to estimate its contribution we can calculate the partition function Z^^^ 
of a single grain (normalized by the partition function of the system without interaction). 
Carrying out the summation over the Matsubara frequencies En and making the gauge 
transformation with the phases (p (r), Eq. (2.15), we obtain 

Z('') = exp/(p), (2.16) 



/(p) = E 



In cosh ~L ^ — In cosh 

2T 2T 



-iNp, 



where is the energy of the state a (the Fermi energy ep is subtracted) . 

The sum over a in Eq. (2.16) extends over all states of the grain. However, the contri- 
bution of the states far from the Fermi energy must be compensated by the local potentials 
N. As usual, an essential contribution comes from energies of the order T in the vicinity of 
the Fermi energy. The linear in p term in the function / (p), Eq. (2.16), must be absent due 
to the choice of N. 

A typical separation of the energy levels is of the order of the mean level spacing 6. 
At not very low temperatures T ':$> 5, the sum over a can be replaced by a proper integral 
over a continuous variable ^ 



E 



and we write the function / (p) in the form 
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1 /-^ 



In cosh ^ — In cosh 

2T 2T 



(2.17) 



where R is an energy in the range T <^ R <^ Sp. Ifp satisfies the inequahty (2.14) we can 
deform the contour of integration when integrating with the first term in Eq. (2.17) and 
integrate first along the straight fine {—R, —R — ip), then along {—R — ip, R — ip) and, at 
last, along {R — ip, R). The integral over the second segment cancels the second term in 
the integrand in Eq. (2.17) and we come to the result 

/W = -^ (2.18) 

Using Eqs. (2.16, 2.18) we can understand easily that the main contribution in integrals 
over p comes from p ~ {TSf^^ < T. Therefore, in the limit T > 5 the variable p can safely 
be put to zero in all calculations and we come to a phase free energy functional containing 
only the phases </> (r) , Eq. (2.15). We see from this argument that introducing the "winding 
numbers" fci, Eqs. (2.15, 2.14), is unavoidable. The opposite hmit T <^ S is more difficult 
and is not considered in the present paper. 

Thus, in the limit T 3> 5, we are able to remove the effective voltage Vi (r) from the 
single grain Hamiltonian Hq, which means removing the Coulomb interaction. However, the 
phases enter now the tunneling Hamiltonian ifj. Expanding in the tunneling term H^, Eq. 
(2.2), up to the second order we obtain the AES^° action S in the standard form 

S^S, + St, (2.19) 

where Sc describes the charging energy 

y 

and St stands for tunneling between the grains 

S, = .9 i: /' drdr-a (r - r') (M^MO^ (2.21, 

|i-j|=a ° \ / 

The function a (r) in Eq. (2.21) has the form 

a (r) = (Re (sin (ttTt + iS))''^)^ . 

In Eqs. (2.19-2.21), i and j stand for coordinates of grains, a is the diameter of a grain. The 
dimensionless conductance g is given by 

g = 27riy^tl (2.22) 

where is the tunneling amplitude from grain i to grain j (spin is included) and is the 
density of states of non-interacting electrons. 
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Calculation of different averages with the functional S, Eqs. (2.19-2.21), implies integra- 
tion over (pi (r) and summation over ki. At large g ^ 1, one can put all ki = 0. However, at 
g % 1 one should sum over all ki and neglecting the contribution of the non-zero winding 
numbers leads to incorrect results. We emphasize that the AES free energy functional can 
be derived under the assumption that the temperatures are not very low, Eqs. (1.3, 1.4), 
and is applicable in this limit only. At lower temperatures one would have to take into 
account interference effects^ or the discreteness of the levels in single grains. 

Using the Kubo formulae, Eqs. (2.5- 2.9), we can express also the conductivity in terms 
of the phases (pi (r) . After simple manipulations we represent the conductivity a {u) as 



10?-'^ 



Jo 



(2.23) 



X(r) = (Xf(r))-^(Xfo(r)XS(0)), 

i 

(r) = e\g J^^ dr' {5 (r) -5{t'- t)) a (r') cos (0i,i+a (r') - 0i,i+a (0)) , 

Xfi (r) ^eng a (r - r') sin (^^i,i+a (r') - 0i,i+a (r)^ dr', 

In Eqs. (2.23), (r) = 0i (r) — 0j (r) for i and j standing for neighboring grains and 

(...) = J (...) exp i-S) D4> y exp {-S) , (2.24) 

where D(p implies both the functional integration over (p (r) and summation over the winding 
numbers ki. 

Eqs. (2.19-2.24) represent the conductivity cr (cu) in a closed form in terms of a functional 
integral. The contribution of the function X2 (r) originates from the diamagnetic term (the 
second term in Eq. (2.6)), whereas the correlation function < XiXi > comes from the 
paramagnetic term 11, Eq. (2.7) (the first term in Eq. (2.6)). 

Although the model described by Eqs. (2.19-2.24) is simpler than the initial model, Eqs. 
(2.1-2.3), explicit formulae can be written only in limiting cases. 

If the temperature T is very high, T ^ Ec ^ ^^^ij^^ where Ec is the electrostatic energy 
of adding one electron to a grain, fluctuations of the phases are negligible and one can 
set = in the expressions for Xi and X2 in Eqs. (2.23). Then, we obtain easily the 
conductivity 



(7o = e^ga^-'^ (2.25) 

Eq. (2.25) describes the classical conductivity of the granular metal without the Coulomb 
interaction and shows that at temperatures exceeding E^. charging effects are not important. 
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In the opposite limit T <C -Ec transport in the granulated system has much more interest- 
ing characteristics. This inequality can be compatible with the inequality T ^ max{(?5, 5}, 
used for the derivation of Eqs. (2.19-2.24), because E^.^ S for 2D and 3D grains (the charg- 
ing energy is inversely proportional to size of the grain, while 5 is inversely proportional 
to its volume). 

At temperatures T < Ec, calculations are possible only in the limiting cases g ^ 1 and 
<^ 1 and this will be done in the subsequent sections. 

The same action S, Eqs. (2.19-2.21), was used in Ref.^^, and a metal-insulator transition 
has been predicted in a 2D array of tunnel junctions. However, the authors of Ref.-*^^ did 
not calculate the conductivity but discussed properties of the partition function. For large g 
they did not account for phase fluctuations properly which, as we show here, are responsible 
for the behavior, Eq. (1.2). Moreover, we find a transition in any dimensionality. 

III. CONDUCTIVITY IN THE METALLIC REGIME AT G » 1 

A. Perturbation theory 

In the limit of large conductances (7 3> 1, the tunnelling term, Eq. (2.21), suppresses 
large fluctuations of 0. It is clear that all non-zero winding numbers ki can be neglected. 
Accounting for non-zero ki (as well as variations of Ni ) would lead to contributions of order 
exp {—g) , which can be neglected in any expansion in 1/g. At the same time, the phase 
fluctuations can change considerably the classical result, Eq. (2.25), even in the limit g ^ I. 

Let us understand first the role of the fluctuations within a perturbation theory in 1/g. 

The zeroth order of the perturbation theory (all phases are set to zero), gives for the 
conductivity the classical result (Tq, Eq. (2.25). In order to consider higher orders we expand 
the action S, Eqs. (2.19-2.21), in 0. 

The quadratic part S2 of the action S will serve as the bare action in the perturbation 
theory we want to develope now. Keeping terms of the second order in in Eqs. (2.19-2.21) 
and performing Fourier transformation in both coordinates of the grains and the imaginary 
time we reduce the action S to the form 

So^TY^ (/)q,„Gq],(/)_q (3.1) 
q,n 

= u;l/ {4E (q)) + 2g ^ (1 - cosqa) , (3.2) 

a 

where E (q) = e^/ (2C (q)) and C (q) is the Fourier-transform of the capacitance matrix Cij 
(q are quasi-momenta for the array of the grains). One should sum in Eq. (3.1) over d unit 
lattice vectors a, where d is the dimensionality of the array. 

If we kept only quadratic in terms in the action S but did not expand the function X2, 
Eq. (2.23), we would reduce the correlation function {X2a {^n)) to the form 

(Xf (Q„))o = ne'g f a (r) (l - e'""^) e-^'^^^^dr, (3.3) 
^0 
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(r) = iTa' J] / -^G,. sin^ f sin^ ^ (3.4) 

where < ... >o means averaging over the phases (j) with the action Sq, Eq. (3.1) and 
Un — 27m. What remains to be done in order to calculate the contribution X2^^ to the 
conductivity is to compute the integral in Eq. (3.3) for the Matsubara frequencies fin and 
make the analytical continuation fi„ —iuo + 5. As it is clear from Eqs. (3.2, 3.4), the 
function Ga (r) contains large logarithms \n{gEcT) and essential r are of the order 1/T. 
Therefore, we may calculate the integral for (t) with a logarithmic accuracy. Neglecting 
the uJ^ term in Eq. (3.2), we reduce Eq. (3.4) to the form 

^ a;„>0 " 

In Eq. (3.5) one should sum over positive Matsubara frequencies up to the cutoff ujc ~ gEc- 
Eq. (3.5) shows a remarkable independence of the result on the structure of the lattice. 
The only information about the lattice is the parameter d entering Eq. (3.5). For the cubic 
lattice considered here, d is the dimensionality of the array. However, for an arbitrary lattice 
the parameter d is equal to the one half of the coordination number. What is also important, 
there are no "infrared" divergencies in the integral over q in any dimensionality including 2D 
and ID. This is specific for the conductivity. We will see later that the tunneUing density 
of states V is sensitive to the dimensionality in the same approximation due to infrared 
divergencies. 

As we are performing the calculations with the logarithmic accuracy, we may replace r 
by 1/T in the function Ga and calculate the remaining integral over r in Eq. (3.3) ignoring 
the dependence of the function G^ on r. Then, we obtain 



< X2 {u) >o= -lue'g j (3.6) 

where 

a = {2T:gdy^ (3.7) 

This form of the correlation function < X2 {ui) > would lead to the power law dependence 
of the conductivity on temperature of the form of Eq. (1.1). A similar result has been 
obtained for the voltage dependence of the conductance of a single junction in a model with 
an electromagnetic environment If Eq. (1.1) were the final result of our calculations we 
could argue that, in the model under consideration, fluctuations of the phases and, hence, of 
the voltages in the grains are equivalent to fluctuations in the electromagnetic environment 
of the works^^'^^. In other words, each grain would be considered as surrounded by an 
effective medium of other grains and the voltage fluctuations of the medium would lead to 
the power law, Eq. (1.1). 

However, Eq. (3.6), is not the flnal result yet because we have to calculate also the 
contribution coming from higher order terms of the expansion of the action S, Eq. (2.19- 
2.21), in (f> as well as the contribution of the correlation function < XiXi >, Eq. (2.23). 
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Taking into account these contributions can be performed writing an expansion for the 
conductivity in powers oil/ g. The first three terms cx*^^^ of the expansion of the conductivity 
in coming from the function < X2 (a;) >o, (3.6) can be written as 

aW/ao = l-«lnf#V-ln^f#l (3-8) 



and we want to find now contributions of the order up to coming from < XiX^ > 
and those originating from 5*4, where 5*4 contains terms of order 0^ in the action S, Eqs. 
(2.19-2.21). 

As concerns a contribution coming from the correlation function < XiXi >, the first 
non- vanishing term is of the order a"^ and it does not contain powers of In (gEc/T). So, 
we neglect the function < XiXi > and concentrate on the contribution to the function 
< X2 (uj) > coming from the anharmonic part 5*4 of the action S. It is clear from a power 
counting that higher order terms of the expansion of S lead to contributions containing 
higher powers of a and we do not consider them now. 

The lowest order contribution coming from S4 is obtained by averaging with the action 
So, Eq. (3.1), of a product of a term 0^ taken from the expansion of X2, Eq. (2.23), and 
S4. A proper expression can be written as 

rl3 rl3 

< V / / / drdndri {e'^"^ - l) a (r) a (n - t[) (3.9) 

X (0i+ao,i (^) - 0i+ao,i (0))^ (0ii+a,ii (n) - 4>h+a,i, {r[)Y >0 

After Fourier transforming the phases in both coordinates and time we can average easily 
with the action 5*0, Eq. (3.1). Then, Eq. (3.9) is reduced to the form 

12T3 Yl Gl,n,G^,,^.\e^^^^° - 1|' n 1^^'''' - 

a,'^>qi>q2,q3 j=l,2 



/ (e^^"^ -1)1 e-'^^'^ -l\^a{T)dT 1 1 - e*'""^^ H 1 - e'^^'^ | (r) (3.10) 
Jo Jo 



X 



Integrals over r can easily be calculated by changing to the variables z = exp (27riTr). The 
lattice integrations are trivial as well and we come to the contribution to the conductivity 
(7^^^ coming from Eq. (3.9) 

„ nin2 2 \ T J 

ni>n2>0 ^ ^ ^ / 

Adding the contributions a^^^ and o"^^^ we obtain for the conductivity a in the limit a; — > 

^ = '^o(l-«ln(^)) (3.12) 
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We see that the terms of the order Xx? {gEc/T) in cr^^^ and cr^^^ cancel each other and the 
accuracy of Eq. (3.12) exceeds aln{gEc/T). We emphasize again that, in this approxima- 
tion, there is no dependence on the structure of the lattice of the grains except that the 
dimensionality (one half of the coordination number) d enters a, Eq. (3.7). However, this 
property holds only for contributions of the type a" In" (gEc/T). Terms with lower powers 
of logarithms depend on the structure of the lattice in a more complicated way. 

The cancellation of the terms of the order In^ (gEc/T) when calculating a, Eq. (3.12), 
is not accidental and we want to demonstrate this within a renormalization group (RG) 
scheme. 



B. Renormalization Group 

In order to sum up the logarithmic corrections to the conductivity we use RG arguments 
suggested for a one-dimensional inverse square XF-model long ago^^ and used later in a 
number of works^°^^^. We assume that the tunnel conductance is large, fif S> 1 and therefore 
we use the phases (f) neglecting the winding numbers. As the starting functional we take the 
tunneling action St 

St^ng T r r drdr'a (r - r') sin^ ( MlljzMll'j (3.13) 
|i-j|=a"^o io V 2 y 

This action contains the conductance g, which determines the conductivity a. The charging 
part Sc is not important for the renormalization group because it determines only the upper 
cutoff of integrations over frequencies. In the limit T — > the function a is proportional to 
(r — t')~^ and the action is dimensionless. 

Following standard RG arguments we want to find how the form of the action St changes 
under changing cutoffs. Generally speaking, it is not guaranteed that after integrating over 
the phases (j) in an interval of frequencies one comes to the same function sin^ (f) in the action. 
The form of the functional may change, which would lead to a functional renormalization 
group. In the present case appearance of terms sin^ 20, sin^ 40, etc. is not excluded and, 
indeed, they are generated in many loop approximations of the RG. Fortunately, the one 
loop approximation is simpler and the renormalization in this order results in a change of 
the effective coupling constant g only. 

To derive the RG equation we represent the phase in the form 

= C + 0ija; (3.14) 

The function (pij^ is not equal to zero in an interval of the frequencies < a; < a;c, while the 
function 0|?2 is finite in the interval Xuc < u < Uc where A is in the interval < A < 1. 
Integrating in the expression for the partition function 

Z = J exp i-St) D(j) 

over the function 0^2, we come to a new action S with the cutoff XcUc- Substituting Eq. 
(3.14) into Eq. (3.13) we expand the action St up to terms quadratic in Integrating 
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over is straightforward and we obtain with the logarithmic accuracy a renormahzed 



effective action St 



where { = — In A. 

We see from Eq. (3.15) that the form of the action is reproduced for any dimensionahty d 
of the lattice of the grains. This allows us to write immediately the following renormalization 
group equation 

^ - - J- (3 16) 

The solution of Eq. (3.16) is simple. Neglecting the Coulomb interaction in the action St, 
Eq. (3.13), is justified only for energies smaller than gE^ and this energy is the upper cutoff. 
Then the renormahzed conductance g (T) takes the form 

and we come to Eq. (3.12) for conductivity. Both the quantities depend on the temperature 
logarithmically. 

Eqs. (3.17) is obtained in the one loop approximation and should be valid so long as 
the effective conductance g (T) remains much larger than 1. Therefore, Eq. (3.12) is also 
more than a result of the perturbation theory and is valid for a/ao ^ l/g. This gives the 
necessary condition for the apphcability of Eqs. (3.12, 3.17) 

1 - {27rdg)-' In {gE,/T) » l/g (3.18) 

However, the condition (3.18) is not sufficient for the applicability of Eqs. (3.12,3.17) 
because the AES action may be used only in the limit 

T>^5 (3.19) 

If the conductance g or the size of the grains are not very large the condition (3.19) can 

become stronger than (3.18). Then, at lower temperatures one should take into account 
interference effects and, depending on the dimensionality d of the array, both metal and 
insulating states are possible. In contrast, Eqs. (3.12,3.17), are valid in any dimensionality. 

Actually, in the main approximation the action St contains only the phase differences 
0ij (r) that can be considered as variables on sites of a lattice dual to the original one. 
These sites are not coupled in this approximation with each other and this explains why we 
obtained the same results as those for one contact. Taking into account the charging energy 
Sc, Eq. (2.20), couples the sites and the structure of the lattice may become important in 
next orders of the renormalization group equations. 
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One can check that the contribution coming from the correlator of the functions Xi in 
Eqs. (2.23) contains additional powers of 1/g and can be neglected in the main approx- 
imation. It is very important that the correlator (XiXi) in Eqs. (2.23) representing the 
"paramagnetic contribution" contains a summation over j, which corresponds to the zero 
quasi-momentum of the function K. Keeping in Xi linear terms in (p would give zero even 
if anharmonic terms in S would have been taken into account. If wc carried out the compu- 
tation for a single grain the contribution from the {XiXi) would not be smaller than that 
coming from {X2). 

It is very important to note that we use the linear response theory for the calculation 
of the conductivity assuming that the external electric field is homogeneous. This contrasts 
calculation of the conductance for a single contact. In our calculations we do not obtain 
a contribution to the conductivity corresponding to the inelastic co-tunnelling known for 
single dots^^. This is natural because in the problem considered the inelastic co-tunnelling 
can occur only through the entire system. This would lead to contributions exponentially 
small in the size of the system. 



IV. TUNNELING DENSITY OF STATES AT G » 1 

The tunneling density of states (DOS) z/; (e) in the grain i can be introduced in a standard 
way through the retarded single particle Green function (e) with both the coordinates 
in the grain i. As we use here the imaginary time representation we calculate first the 
temperature Green function (r) at Matsubara frequencies = nT {2n + 1) . This leads 
to a function i>i (£„) 

z>i (£„) = -TT-i J dre^'-^g,, (r) (4.1) 

The DOS Vi [e] can be found by the analytical continuation 

i/i(£) = Im(z>i(£„)|,^_,,J (4.2) 

Following the same procedure as the one used in the previous sections we perform the gauge 
transformation, Eq. (2.11), reducing the calculation to integration over the phases 0i(r). 
As a result, we obtain for the function z/j (£„) 

V, (£„) = v^T j^^ '^^^^ (^""P ~ ) ^^-^^ 

where the symbol (...) means as before averaging with the action Eqs. (2.19-2.21) and 
integration over 0; (r) includes summation over the winding numbers ki. The principal value 
of the integral over r is implied in Eq. (4.3). 

Galculation of DOS is quite different from that for the conductivity because Eq. (4.3) 
contains the phases 0i (r) corresponding to the grain i but not the phase differences 0ij (r). 
This can lead to a non-trivial dependence of the final result on the dimensionality d of the 
lattice. 
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In the limit of large 5' ^ 1 wc expand, as when calculating the conductivity, the functional 
Sj Eqs. (2.19-2.21), in 0ij (r). If the d = 3 one should expand in (r) both the action S and 
the exponential in Eq. (4.3). The reason is the same as for calculation of the conductivity: 
all terms of the expansion give additional logarithms and there is no reason to keep in the 
action only quadratic terms and at the same time not to expand the exponential in Eq. (4.3). 
Therefore in the 3D case one can expect logarithmic corrections to the tunneling density 
of states with coefficients depending on the structure of the lattice. Making analytical 
continuation onto real energies s we can write the tunneling density of states (s) in 3D 
in the main approximation as 

where 



U!n>0 ' 



and Gqn is given by Eq. (3.2). 

Calculating the sum over cun with the logarithmic accuracy and making the analytical 
continuation we obtain 

(s) = f 1 _ A In f^l^l—) ) (4.6) 

where 

^_ f a^d^q 1 ,^ „x 

~ J {2nf Ea(l-cosqa) ^ ' 

Eq. (4.6) shows that the density of states (e) of the 3D array of the grains has a 
logarithmic dependence on temperature like the conductivity cr, Eq. (3.12). However, in 
contrast to the latter, Eq. (4.7) is valid only if the logarithmic term is much smaller than 1. 

The situation is more interesting in one and two dimensional systems. In this cases the 
integral over q in Eq. (4.7) formally diverges, which means that one should take into account 
the u;^ term in the function Gq„, Eq. (3.2). This term cuts the infrared divergency in the 
integral over q but one obtains in the density of states a stronger singularity in e, T than 
the one in 3D, Eq. (4.6). 

Fortunately, this makes the calculation even easier and allows us to obtain explicit results 
in the non-perturbative regime when the DOS considerably deviates from uq. This simplifi- 
cation is due to the fact that the strongest singularities come from the expansion in 0i (r) 
of the exponential in Eq. (4.3). The anomalous contributions originating from small q arise 
in the expansion of 0i but cancel each other in the contributions coming from expansions 
in 0ij (r) = 0i (r) — 0j (r) in the action Sf, Eq. (2.21). Therefore, expanding in 0ij (r) one 
obtains usual logarithmic contributions of the type {gin {gEc/T))^ , which allows us to keep 
in the action S only quadratic in 0ij (r) terms. 

Then the integration over (r) in Eq. (4.3) can easily be carried out and we obtain for 
the low dimensions 
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r0 lEnT 

^ (^) = ^oT ^^^^ exp (-G (r)) (4.8) 

It is clear that the singularity in the exponent in Eq. (4.8) is stronger than logarithmic and 
this justifies the approximation used in the derivation. 

In the most interesting case of a 2D array of the granules both the summation over the 
frequencies and integration over the momenta q give logarithms and we come to the final 
result 

-(^)-oexp(-^ln^(;jjf|^)) (4.9) 

Eq. (4.9) is valid down to max(£,T) ~ g5 when the description in terms of the phase 
functional Eqs. (2.19-2.21), is still applicable. Eq. (4.9) perfectly agrees with the 
corresponding result obtained long ago for disordered films using a replica cr-model^. This 
result was reproduced for disordered films in a number of subsequent publications^^"^^ using 
different approaches. The strong anomaly in the exponent in Eq. (4.9) is due to the fact 
that the one particle Green function is not gauge invariant^. The singularity is formed by 
almost pure gauge fluctuations of the electric fields. Gauge-invariant characteristics like 
conductivity are not influenced by such fluctuations and therefore are less anomalous. 

The applicability of Eq. (4.9) not only for disordered fllms but also for the granular 
systems at not very low temperatures shows that the result is very robust. In contrast, 
the dependence of the conductivity on temperature, Eqs. (3.12-3.17), cannot be used for 
very low temperatures or disordered "homogeneous systems". The formal reason for this 
difference is quite clear: the main contribution to the conductivity comes from momenta 
q ~ a^^, whereas the main contribution to the density of states comes from small q <^ a^^ . 
The latter hmit is not sensitive to the structure of the system at short distances. 

The coefficient in front of In^ is somewhat different than that of the recent work^^. This 
is because we assumed that E (q) remains flnite in the limit of g — > (there is a screening 
in the system). In the integral over q in Eq. (4.5, 3.2), essential {qaf were of the order of 
uj/ gE^. In contrast, if one starts with a non-screened 2D Coulomb interaction Vq = 27re^/g, 
the essential q are proportional to cu and this increases the coefficient by the factor of 2. 

The first term of the expansion of the exponential in Eq. (4.9) is just the Altshuler- 
Aronov correction generalized to the case of the granular metal^. Of course, the same is true 
for the correction to the conductivity, Eq. (3.12). This is quite natural because at not very 
low temperatures weak localization corrections are suppressed but the Altshuler-Aronov 
corrections still give important contributions. Actually, the function Icu^ICq, Eq.(3.2), is 
just the Goulomb propagator screened by the electron-electron interaction. In other words, 
the theory developed in the previous sections for g ^ 1 starting from the AES action is 
another way of calculation of the Altshuler-Aronov corrections for the granular systems. 
These calculations could be performed diagrammatically, although establishing the non- 
perturbative results for both the conductivity, Eq. (3.17), and the DOS, Eq. (4.9), would 
be considerably more difficult. Clearly, the action S described by Eqs. (2.19-2.21) may 
not be used at zero temperature and the dissipation resulting from this action is not a zero 
temperature effect. 
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If the tunnel conductance g becomes of the order of 1 or smaller diagrammatic expansions 
are no longer helpful and calculations without the phase functional, Eq. (2.19-2.21), are 
hardly possible. In this regime a proper account of non-zero winding number ki is very 
important. In the next section we show how calculations can be carried out in the limit of 
weak coupling constants g' <^ 1. 



Study of the granular system described by the action S, Eqs. (2.19-2.23), for an arbitrary 
conductance g is difficult. The logarithmic behavior, Eq. (3.17), describes the conductivity 
of the granular system at sufficiently large (? > 1. Eqs. (4.6-4.9) are also applicable only in 
this limit. Not being able to consider the model for arbitrary g, we restrict ourselves with 
the limit of small g <^ 1. 

We will see that the temperature dependence of both the conductivity and the DOS 
becomes exponential in this limit. This means that increasing the tunneling amplitude 
at a fixed T we go from the almost metallic regime to an insulating one. This can be 
achieved experimentally changing the coupling between the grains while measuring at the 
same temperature^. 

Calculation of physical quantities at small g <^ 1 can be performed expanding the func- 
tional integral in Eq. (2.24) in the tunneling part St, Eq. (2.21), of the action. When 

calculating conductivity the main contribution comes again from the function {X2 (r)) in 
Eqs. (2.23). As concerns the tunneling density of states we can use as before Eq. (4.3). In 
the lowest order one can completely neglect St in the both formulae. Then, the calculation 
of the DOS, Eq. (4.3), reduces to computation of the phase correlation function 11 (r) 



where the phases 0i (r) are introduced in Eq. (2.15), Sc is given by Eq. (2.20) and the 
averaging should be performed with this functional. As concerns the conductivity, we need 
a slightly different correlation function n (r) 



The phase correlation function 11 (r) is somewhat simpler and let us show in detail how to 
calculate it. A proper modification for 11 (r) is simple. 

The computation of the average in Eq. (5.1) can be performed using two different 
methods. A more straightforward way of calculating is to use the definition of (pi (r), Eq. 
(2.15), which allows us to represent the action Sc as 



V. WEAK COUPLING BETWEEN THE GRAINS 



A. Phase correlation function 




(5.1) 




(5.2) 




(5.3) 




(5.4) 
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S,[k]^T7r'Y,h{B-')^.k^ (5.5) 

ij 

where 

Writing Eqs. (5.4-5.5) we neglected integration over the variables pi from Eq. (2.14). As we 
have discussed in Sec. II, in the limit T ^ S the main contribution comes from p ~ (T5)^^^ 
and we can simply put in all expressions p — 0. Using Eqs. (5.3-5.5) one can carry out 
integration over the phase and summation over the winding numbers separately. The 
phase correlation function 11 (r) can be written as 

n (r) = (exp(-i (0; (r) - 0; (0)))^ (exp(-27riA;iTr)), (5.6) 

Integrating over the phase (r) we obtain for < r < /3. 

(exp(-z (0i (r) - 0i (0)))^ = exp (-Bn (r - Tr^)) , (5.7) 

Eq. (5.7) was used for the function n (r) in many previous works^^'^^. However, Eq. (5.7) 
is not the final result because the function 11 (r), Eq. (5.6), contains the second average and, 
in addition to integrating over 0i (r) , one must sum over the winding numbers ki. 

Calculation of the second average in Eq. (5.6) can be performed using the Poisson 
formula 

k=—oo m=—oo °° 

for any function / (x) . 

As a result, we rewrite the second average in Eq. (5.6) as 

(exp(-27ri/.iTr)), = E / ( "J E + ^ E (^J " ^^'^'j) ) 

{mi}-^ V jl j / 

(5.8) 

where Z is the normalization factor that can be obtained putting r = in the integral in the 
second line of Eq. (5.8) and T>x = Yl^dxi. Summation over all integers rrii for all grains is 
implied. The integration over Xi should be performed in the infinite limits and the integral 
in Eq. (5.8) can be easily calculated. Substituting the result of the integration into Eq. 
(5.6) and using Eq. (5.7) we find for the phase correlation function 11 (r) 

n (r) = Z-^ exp {-But) E ( " E 2^^k5ki - /3 E ^kimtmi j (5.9) 

{mj} V k k,l / 

where all mj are integers and Z is the normalization coefficient (11 (0) = 1). The quantity Ba 
in Eq. (5.9) is the charging energy of an extra electron in a grain i in an otherwise neutral 
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system without excitations. The necessary periodicity in r of the function IT (r) with the 
period /? is evident from Eq.(5.9). The final resuh for the phase correlation function 11 (r), 
Eq. (5.9), is essentially different from Eq. (5.7) obtained by neglecting the contribution of 
non-zero winding numbers. The exponent contains only linear in r terms, which contrasts 
Eq. (5.7). 

Actually, the form of Eq. (5.9) is absolutely natural and can be obtained using the 
standard quantum mechanical formalism instead of calculating the functional integrals in 
Eq. (5.1). Within this formalism velocities in the action in the functional integral should 
be replaced by the corresponding momentum operators in the Hamiltonian. So, instead 
of having the derivatives d(f)/dT in the action Sc, Eq. (2.20) one would have to write the 
operators d/d(f). The presence of winding numbers in the action introduces periodicity. As a 
result, the corresponding angle variables in the Hamiltonian formulation should be taken 
in the interval [0, 27r] and all wave functions of the Hamiltonian must be periodic in 4> with 
the period 27r. 

As the result the calculation of the phase correlation function H (r), Eq. (5.1), reduces 
to calculation of quantum mechanical averages with the effective Hamiltonian Hes 

Hes = Bijpipj, pi = -id/d(pi, (5.10) 

ij 

where the "angles" <pi vary between and 27r. The correlation function H (r) can be written 
in the form 

H (r) = (e-<^'(")-^'(°)))^^^, 0i (r) = e^««>ie-^««^ (5.11) 

In order to calculate the average with H^^ in Eq. (5.11) one should find eigenfunctions 
^{„}0f the Hamiltonian ^eff, Eq. (5.10). These eigenfunctions have the simple form 

'^{n}^\{eK^{ini(t)i) (5.12) 

i 

where n\ are integers. Calculating the matrix elements entering Eq. (5.11) and performing 
summation over all states with the weight exp (— i^i^j/T), where -E{n} are eigenenergies of 
the Hamiltonian ^eff, Eq. (5.10), we come to Eq. (5.9). 

The operators pi and 0i are conjugate to each other with the commutator 

[ph 0i] = -i 

Therefore we can alternatively write the operators 0i as 0i = d/dpi. The eigenvalues of the 
operator pi are integers. The operator exp ^±z<^j acts as follows 

exp {±4) f (p) = exp i^d/dp) f{p)^f{p± 1) (5.13) 

This gives another convenient way of calculation of the quantum mechanical averages. 

The use of the Hamiltonian formalism for calculation of functional integrals over the 
phase (f) of the superconducting order parameter has been suggested in an earlier work 
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on granulated superconductors^^ . Within this approach the effective Hamiltonian Hes, Eq. 
(5.10), was derived for the operator p of the number of Cooper pairs and the phase correlation 
function IT (r), Eq. (5.9) has been obtained. As concerns the normal metals the correct form 
of the phase correlation function 11 (r), Eq. (5.9) has been written for the first time in our 
previous paper^^. 

The present consideration demonstrates explicitly that accounting for the winding num- 
bers leads to the charge quantization. The function II{t), Eq. (5.2) can be calculated in 
the same way and can be written as 

n (r) = exp i-B^r) exp | - 2rmk (5k,i+a - ^ki) - ^Ji^J^i ) (5-14) 

where = -Bii + i?i+a,i+a — 5i,i+a — -Bi+a,i is the energy of an electron-hole excitation between 
neighboring grains, a = j — i. As the array of the grains is regular we omit the subscripts 
i, j when writing the correlation functions 11 (r) and 11 (r). 

The functions 11 (r) and IT (r) allow us to calculate the tunnehng density of states and 
the conductivity using Eqs. (4.3) and (2.23). 



B. Tunnelling density of states at g ^ 1 

Substituting Eq. (5.9) into Eq. (4.3) we have to calculate the remaining integral over r 
and perform the analytical continuation — > —ie + S. It is convenient to shift the contour 
of the integration and integrate first along the line (0,ioo), then along {ioo,ioo + (5) , and 
finally along (ioo -|- /3, /?). The integral over the second segment vanishes and we reduce the 
function Di (£„), Eq. (4.3), to the form 

V, {e) = - 2T r dt^^^ilM sin (t5„) (5.15) 
Jo smh {nTt) 

xZ~^ exp(— /? mkBy^irrii — 2it m^Bki)] 

{mi} kl k 

Now the analytical continuation can be done easily. Taking the imaginary part of the 
function u according to Eq. (4.2) and using the integral 



sm au 1 a 

du — - tanh - 



/o sinhTTU 2 2 

we obtain for the density of states u (e) 



Z-^exp(- g^'"'j;^-'"" ) (5.16) 



° {"Ik} 



where 
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Z = ^ exp ( _ ^ki--K-K..-. J ^5^^^ 



{"^k} 
and 



n[x) — 



+ 1 

is the Fermi distribution function. In Eqs. (5.16, 5.17) summation over all positive and 
negative integer rrik should be performed. 

The function u (e) is even in the energy e and approaches 1 in the limits T — > oo or 
l^l oo. At low temperatures T ^ B^, \6\ and \e\ < Ba the main contribution comes from 
the ground state configuration when all mi = 0. In this limit we obtain 

u(s) f B\ F , , 

2exp -- cosh- (5.18) 



where B = Ba. Eq. (5.18) demonstrates that there is a gap in the density of states and 
this gap is equal to the single electron charging energy of the grain. Of course, Eq. (5.16) 
corresponds to a fixed chemical potential in the grain, which means that the grain is not 
completely isolated and there are processes that keep the chemical potential fixed. 

One can easily generalize Eq. (5.16) to the case when an additional voltage is present in 
each grain. Writing the effective chemical potential /li of a grain i as 

j 

we can generalize Eq. (5.16) by replacing everywhere mk + -/Vk- If -^k is randomly 

distributed over the grains one should carry out an additional averaging over this variable. 
This is beyond the scope of the present paper. 

To the best of our knowledge, Eq. (5.16) has not been written before, although it clearly 
follows from the "orthodox theory" of the Coulomb blockade^^ . The latter can be seen from 
the fact that the one particle part Hq of the Hamiltonian, Eq. (2.1) commutes with the term 
He, Eq. (2.3), describing the charging energy. This means that in the absence of tunneling 
between the grains the charging energy should be simply added to the Fermi energy for 
non- interacting particles. Eq. (5.16) clearly corresponds to this picture. For an arbitrary 
configuration {mt} of charges on the grains the energy of adding one electron (hole) to a 
grain is equal to 

5±2^mk5ki 

k 

We see that just this energy enters the Fermi distribution functions in Eq. (5.16) shifting 
the Fermi energy. The exponential in Eq. (5.16) is the weight for the configuration of the 
charges in the system and two Fermi functions describe contributions of electrons and holes. 
Although the direct derivation "in the electron language" would be simpler, the present 
calculation demonstrates explicitly how the phase functional Sc, Eq. (2.20), works. We 
emphasize that without the summation over the winding numbers the correct result could 
not be (and has not been) obtained. 



22 



C. Conductivity at g 



Calculation for the conductivity can be performed in the same way as for the tunnelling 
density of states. It is important that in the limit g <C Ithe main contribution comes from 
the term X2 in Eq. (2.23). The contribution coming from the term Xi is of higher order in 
g and we neglect it. Then, the Fourier transformed response function K (Qn), Eq. (2.23), 
takes the form 

K{n^) = 27re'gT' [' dr^-^^^P^flir) (5.19) 
Jo sm^ ttTt 

where, again, the principal value of the integral should be taken. 

Shifting the contour of integration over r in the same way as it has been done when 
calculating the tunnelling density of states we can perform the analytical continuation fl^i 
—iuo + 5. In the limit of small frequencies uj we expand the integral in uj. The first non- 
vanishing term is linear in uj and we reduce the formulae for the conductivity a using Eqs. 
(2.23, 2.25) to the form 

2nT dt—^^^^ sin {tBf) Cf {t)^ (5.20) 

where 

Cf {t) = Z-^ ( -(^ H ^kSkimi j cos ( 2t J] mk (Sk,i+a - Ski) ] , (5.21) 

{mk} \ kl / \ k / 

Z is given by Eq. (5.17) and gq is the classical conductivity, Eq. (2.25). Calculation of the 
integral over t in Eq. (5.20) can easily be performed using the formula 

[°° , usiniau) 1 ,^ ^ / NX 

/ • .2 , (5-22) 
Jo smh (ttw) 27r 



/(a)=2 "--'>--°+r 
(1 — e ") 



and we obtain finally 



C7 = (Jo Z-' J2 "^k^kim, / f3{Bt + 2 5^ mk {B^,i+^ - Ski)) (5.23) 

{mk} V kl / \ k / 

We see from Eq. (5.23) that, in order to get an explicit expression for the conductivity, 
one should sum again over all charge configurations (In Eq. (5.23) the conductivity is 
calculated between the grains i and i + a). 

The limit of high temperatures T exceeding the charging energies B can be obtained using 
the property /(O) — > 1 when a — > 0. In this limit we come to cr = a"o, which demonstrates 
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that as soon as the Coulomb energy is not important the transport is described by the Drude 
formula. 

In the opposite limit, T -C i?, the main contribution in Eq. (5.23) comes from charge 
configurations with the lowest charging energies. For the calculation of the conductivity we 
need also the asymptotics of the function / (a) in the limit a — > oo. In this limit we may 
write this function as 



As when calculating the tunnelling density of states, we consider first the contribution 
of the ground state configuration with all mi = 0. However, in contrast to the tunnelling 
density of states, this configuration gives not necessarily the main contribution because it 
contains a large energy of a dipole consisting of an additional electron in the grain i 
and a hole in the grain i + a (or vice versa). The energy of this dipole is equal to 2B if we 
neglect non-diagonal Bij. If non-diagonal components Bij are not equal to zero the dipole 
energy is smaller but we assume that it is larger than the charging energy of one electron 
B. Physically, this contribution corresponds to transport in a completely neutral grains. 
In order to contribute to the current an electron must jump from one grain to another. 
However, this costs an energy which is just the dipole energy 5?. So, the contribution of 
the ground state configuration can be estimated as exp {—Bf/T) and we want to show that 
a larger contribution exists. 

The most efficient process contributing to the current is when an additional charge exists 
in a grain but all other grains are neutral. Then, jumping from grain to grain costs no 
energy. As such a configuration is not the ground state, the probability to have this state 
is proportional to exp {—B/T). However, the overall contribution to the conductivity is in 
this case larger because B < Bf. 

This picture clearly follows from Eq. (5.23). Two configurations give the main contribu- 
tion to a in Eq. (5.23). We can put = 1 at k = i and = for all k 7^ i or = — 1 
for k = i + a and mk = for all k 7^ i + a. In both the cases the argument of the function 
/ in Eq. (5.23) is equal to zero and we obtain for the conductivity at low temperatures 



Eq. (5.25) shows that in the limit of small couplings between the grains g the macroscopic 
conductivity is exponentially small in temperature. This is a typical example of an activation 
process. 

Of course, the exponential behavior of the physical quantities, Eqs. (5.18, 5.25), was 
derived under the assumption that all the grains are mesoscopically equal (they have the 
same size and shape but may have small irregularities different for different grains). In real 
samples the shape and the size may vary and qualitative estimates show that instead of the 

activation law a dependence of the type exp I —A/-\/t\ can be more proper for this case^^. 



We studied effects of the Coulomb interaction on the conductivity a and tunneling density 
of states v (e) of granular metals. Calculations with the Hamiltonian for interacting electrons 



/ (a) ~ 2a exp (—a) 



(5.24) 



(7 = 2(7oexp(-B/r) . 



(5.25) 



VI. DISCUSSION 
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were reduced to calculation of functional integrals with a phase action of a form proposed by 
Ambegaokar, Eckern Schon^°. This action has been derived recently for the granular systems 
microscopically^, which allowed to clarify conditions for its applicability. These conditions, 
Eqs. (1.3, 1.4), correspond to the limit of not very low temperatures, such that weak 
localization effects are suppressed. In the limit of large tunnelling conductances g the results 
obtained with AES functional correspond to Altshuler-Aronov corrections^ and could be 
calculated diagrammatically (although this way of calculations would be more complicated). 
At smaller g non-zero winding numbers become very important and we developed a proper 
scheme of calculations. 

Although the interference effects leading to locahzation corrections are neglected at such 
temperatures, interesting effects occur. In the limit of large g a logarithmic dependence of the 
conductivity a on temperature, Eq. (3.12), is obtained. Eqs. (3.12, 3.17) are applicable in 
any dimension and at any magnetic field, which distinguishes it from the weak localization 
correction. The logarithm in Eqs. (3.12, 3.17) is not just a small correction and these 
formulae are applicable until the conductance becomes of order unity (It is important that 
the conductivity and not the resistivity is linear in the logarithm of the temperature). 

In contrast to the conductivity a, the tunnelling density of states u (e) is dependent on 
the dimensionality d of the system in a non-trivial way in the limit g ^ 1. In 3D it has a 
logarithmic correction, Eq. (4.6), but is described by a more complicated formula, Eq. (4.9), 
in 2D. This formula is well known for disordered films^'^^ where it was obtained within the 
(T-model approach. We see that the same formula is valid for the granular metal at not very 
low temperatures and it can be obtained without cr-models. 

In the limit of a low coupling between the grains both the tunnelling density of states 
and the conductivity are exponentially small, Eqs. (5.18, 5.25). This is due to a finite 
charging energy arising when the electron tunnels from one to another grain. Although the 
results in the limit of vanishing g are rather simple, it was important to derive them from 
the AES phase functional with a proper summation over the winding numbers. This way 
of calculations is very close to the one suggested previously for granular superconductors'^^. 
In some experiments^ the dependence of the conductivity on temperature is described not 
by the activation type formulae but by the function exp (^—a/^/T^ . This dependence may 

originate from fluctuations of the charging energy of the grains^^. 

Comparing Eqs. (3.12, 4.6, 4.9) obtained in the limit g :$> 1 with Eqs. (5.18, 5.25) derived 
for g '^1 we conclude that there must be a considerable change of the temperature behavior 
of the physical quantities when changing the coupling between the grains. Of course, this 
cannot be a sharp transition because we consider the limit of finite temperatures and cannot 
extrapolate the results to T = 0. However, tuning the coupling experimentally at a given 
temperature one may see the change of the regimes that would look like a "metal-insulator" 
transition. We emphasize that this "transition" should occur in any dimension of the array 
of the grains. 

It is not clear from the present consideration whether there should be a sharp transition 
from the metallic to the insulator state at a critical value gc and this is a definitely interesting 
problem for a further investigation. This problem is closely related to the question whether 
the activation energy (Coulomb gap) turns to zero or has a jump at g = gc- 

The model of the granular metal may be used to describe disordered electron systems 
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at low electron density. In such systems electrons can spend a considerable time in traps 
or "puddles" that can be due to strong fluctuations of a disorder potential. Considering 
such systems with the model of a granular metal may be a reasonable approximation. In 
this case, potential wells where electrons are trapped would correspond to the grains in the 
model of the granular metal. 

The model considered in the present paper and the results obtained can be relevant to 
many experiments on different materials. Of course, specially prepared granular metals like 
those considered in Refs.^'^^'^^ should be the first object of the application of the theory 
developed. As the model and the results obtained are quite robust, one can expect that the 
corresponding phenomena have been observed in the granular materials. 

This is really so and many experimental data can be explained in the framework of our 
model. First, let us make a comparison of our theory with experimental results of Ref.^ 
on films made of Al grains embedded in an amorphous Ge matrix. At low temperature 
superconductivity of in Al grains was destroyed by a strong magnetic field. Depending on 
the coupling between the grains (extracted from the conductivity at room temperatures) 
the samples of the experiment'' were macroscopically either in an insulating state with the 
temperature dependence of the resistivity R ~ exp (a/T^/^) or in a "metallic" one. However, 
the resistivity of the metallic state depended on temperature and the authors suggested Eq. 
(1.1) to describe this dependence. As the exponent a for the "metallic" sample was small 
we may argue that Eq. (1.2) should not be worse for fitting the experimental data. Then, 
we can estimate the exponent a without using fitting parameters. 

The sample of the experiment'' had the room temperature resistivity Rq = 7.3 x 10^^r2cm. 
The diameter of the grains was 120±20A, which allows, using the value h/e^ = 4.1 x lO^Q, 
to estimate the dimensionless tunnel conductivity as g = 0.7. If we put d = 2 in Eq. (3.7) 
we obtain a — 0.116, which would perfectly agree with the experimental value a — 0.117 
from Eq. (1.1). However, everything is not so simple because the films used in Ref.^ were 
rather thick and, at first glance, one should use d = 3. This would change the result by 30% 
making the agreement less exciting. Nevertheless, the value of d in Eq. (3.7) corresponds 
rather to the half of the contacts of a single grain than to the real dimensionality. Therefore 
the experimental value of a indicates that either the grains are not closely packed such that 
the typical number of contacts per grain is 4 or our calculation is too rough to provide a 
quantitative agreement with the experiment (the value of a, Eq. (3.7), is based on the 
assumption g ^ I but the experimental value of g is of order 1). 

The resistivity of samples with a high room temperature resistivity (a weak coupling 
between the grains) behaved as exp (a/T^/^) rather than obeying the activation law, Eq. 
(5.25). According to Ref.^^ this can be attributed to a variation of the size of the grains or 
of the local potential. However, a model considered in Ref.^^ is rather special because the 
distance between the grains and the tunnelling amplitude was assumed to be related to the 
charging energy in a certain way. The law exp (a/T^/^) is rather common for granular mate- 
rials with weak couphng (see also e.g. Ref.^^) and the reason for such an universahty is still 
not clear. A random hopping conduction mechanism of Ref.^^ suggested for semiconductors 
can hardly be used for the granular metals. 

A logarithmic dependence of the resistivity on temperature has been observed in other 
granular materials. In Ref.^^ a granular cermet consisting of NbN grains in a boron nitride 
insulating matrix was studied. Again, at small coupling between the grains the temperature 
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dependence of the resistivity exp (a/T^/^) was observed in a very broad interval of temper- 
atures. The resistivity of samples with a strong coupling between the grains was very well 
described by the law 

R^Roln{To/T) (6.1) 

which is close to Eqs. (1.1, 1.2) if the temperature interval is not very large such that the 
variation of the resistivity is small. However, the law, Eq. (6.1), gave a good description for 
the temperature dependence of the resistivity in a very broad region and the changing of the 
resistivity was not small. The reason for the applicabihty of Eq. (6.1) in a so broad interval 
of temperatures is not clear because according to the results of the renormalization group 
analysis of Sec. Ill not the resistivity but the conductivity should obey Eq. (6.1). A more 
careful experimental study might clarify this question. Anyway, the logarithmic behavior of 
Refs.^''^^ remained unexplained at all and our work is the first attempt to construct a theory 
of this effect (an explanation in terms of weak localization corrections or the Kondo effect 
can be excluded immediately because the logarithmic temperature dependence was observed 
also in very strong magnetic fields and the systems were three dimensional). 

The unusual logarithmic behavior of the type, Eq. (6.1), has been observed not only 
in "standard" granular systems but also in high-Tc cuprates at very strong magnetic fields. 
The first observation of this dependence was done on underdoped La2-xSrxCuOi crystals^^. 
The superconductivity in this experiment was suppressed with pulsed magnetic fields of 61T. 
It was found that both the in-plane resistivity pah and out-of-plane resistivity pc diverged 
logarithmically with decreasing the temperature. This means again that a 3D effect was 
observed in a very strong magnetic field and traditional explanations like localization or 
Kondo effect could not clarify the situation. 

In a subsequent publication^^ a metal-insulator crossover was observed in the same mate- 
rial at a Sr concentration near optimum doping [x ~ 0.16). In underdoped samples both pab 
and Pc showed no evidence of saturation at low temperatures and diverged as the logarithm 
of the temperature. The authors called this state "insulator" in contrast to the state at high 
doping where the resistivity did not have a pronounced dependence on the temperature. It 
was conjectured in Ref.^^ that the logarithmic behavior they observed might be related to 
the one seen in the experiment^'^ on granular NhN . This would demand a phase segregation 
throughout the underdoped regime of LSCO. However, as no explanation had been given 
for the logarithmic behavior in Ref.^^, no explanation has been given to the experiments^^'^^ 
either. We hope that our results for the model of the granular materials may be appli- 
cable to the experiments on the La2-xSrxCu04 crystals^^'^^, which would mean that the 
underdoped crystals have a granular structure and the logarithmic behavior is due to the 
Coulomb interaction. The transition to the metallic state of Refs.^^'^^ would mean that at 
higher doping the granularity disappears. 

The logarithmic dependence of the resistivity on temeperature has also been observed in 
many other experiments. For example, in Ref.^^ this dependence was observed in granular 
Pb films. It was also observed in phase compounds of Nd2-xCexCu04^y, Ref.^^. In each 
case the reason for such a behavior was not clear. 

An interesting conclusion has been made recently about the structure of a 2D gas in 
GaAs/AlGaAs heterostructures where a metal- insulator transition was observed^^. Mea- 
suring the local electronic compressibility the authors found that the metallic phase was 
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homogeneous in space, which is naturaL In contrast, the system becomes spatially inhomo- 
geneous as it crosses into the insulating phase. The structure seen in the insulating state 
indicates that the electrons were located in "puddles". So, modelling the system in terms 
of a granular metal might be a reasonable approximation also in this case. 

We are grateful to A. Altland and B.L. Altshuler, A.V. Andreev, L.I. Glazman, A. 
Gerber, A.I. Larkin, and Yu.V. Nazarov for discussions. A support of the SFB/Transregio 
6029 and German-Israeli Foundations (GIF and DIP) is greatly appreciated. 
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